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Abstract 

As  a  two  year  project  related  to  our  earlier  AFOSR  work,  we  examined  the  use  of  surface 
textures  to  suppress  the  growth  of  turbulent  spots  during  the  later  stages  of  boundary 
layer  transtion.  The  textures  are  small  and  closely  spaced  and  hence  require  detailed  direct 
numerical  simulation  of  the  near  surface  flow  to  capture  the  physics.  The  work  emphasized 
the  fundamental  nature  of  that  fiow/texture  interaction.  This  project  primarily  involved  the 
application  of  our  pre-existing  computational  techniques  and  software  but  did  require  the 
alteration  of  the  code  to  handle  a  boundary  layer  and  to  run  with  high  spatial  resolution 
on  a  parallel  processor.  The  possible  applications  to  boundary  layer  control  are  suggested 
iu  Background  while  the  section  called  Proposal  Objectives  reviews  the  approaches 
considered.  Finally,  an  attached  Detailed  Report  is  a  comprehensive  discussion  of  previous 
scientific  work,  our  techniques,  and  our  results. 

Background: 

Boundary  Layer  Transition:  A  greater  knowledge  of  turbulent  boundary  layer  structure 
and  development  and  the  creation  of  means  for  boundary  layer  transition  control  are  essential 
to  several  applications  in  the  fields  of  aerodynamics.  As  Tollmien-Schlichting  waves  in  a 
laminar  boundary  layer  grow  to  develop  unstable  3D  perturbations,  turbulent  spots  form.  In 
an  uncontrolled  flat  plate  boundary  layer,  such  arrowhead- shaped  spots  appear  at  randomly 
distributed  points  and  their  growth  and  merging  lead  to  complete  transition  to  turbulence. 
A  spot  appears  as  a  fairly  well  defined  region  of  turbulence  surrounded  by  laminar  flow. 
Wygnanski  e.t.  al  [1].  performed  detailed  experimental  studies  of  such  spots  and  Leonard  [2] 
briefly  reviewed  some  of  the  first  simulations  with  vortex  methods  of  a  spot- like  perturbation 
iu  a  laminar  flow.  Spots  have  been  sucessfully  simulated  in  several  shear  flows  (e.g.,  [3]). 
Spots  appear  to  displace  some  near-wall  fluid  outward  (a  spamvise  velocity  component)  and 
produce  spamvise  traveling  waves  [4]. 
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Often,  flow  conditions  are  not  sufficiently  perfect  that  transition  occurs  in  the  classic 
2D  TS  instability  — *  3D  perturbation  growth  — *  turbulent  spot  progression.  For  example, 
roughness  elements  (e,g.,  bugs)  can  cause  a  premature  immediate  bypass  transition  to  tur¬ 
bulence  in  a  wedge- like  region  aft  of  the  perturbation.  (Although  this  situation  might  be 
conjectured  to  be  just  a  continuous  sequence  of  spots.)  Joslin  [5]  particularly  points  out 
the  detrimental  effects  of  leading  edge  contamination  for  laminar  flow  control.  A  turbulent 
boundary  layer  is  characterized  by  swirling,  unsteady,  seemingly  random  motion.  In  fact, 
however,  the  motion  is  not  random  but  may  be  organized  into  a  variety  of  quasi- periodic 
vortical  structures.  These  vortices  lead  to  the  formation  of  low  speed  streaks  which  are  be¬ 
lieved  to  be  responsible  for  the  breakdown  and  eruption  of  the  wall  layer,  a  key  mechanism 
in  the  production  of  new  vortical  structures  and  the  increase  in  shear  stress  at  the  wall. 
We  note  that  in  both  the  formation  of  turbulent  spots  and  wedges,  the  near-wall  turbulent 
structures  of  streamwise  vortices  and  hairpins  appear  much  like  those  in  a  fully  developed 
turbulent  boundary  layer.  In  fact,  the  lateral  spread  with  downstream  propagation  distance 
of  those  spots/wedges  appears  related  to  t  hose  same  fundamental  structures. 

If  one  can  delay  or  confine  the  regions  of  laminar  to  turbulent  transition,  the  benefits 
could  be  appreciable.  For  example,  since  the  viscous  drag  of  a  turbulent  boundary  layer  can 
be  four  or  more  times  greater  than  that  of  a  laminar  layer,  one  can  imagine  that  sufficiently 
delaying  transition  could  lead  to  greater  aircraft  ioiter  times  over  targets  or  a  decrease  in  the 
need  for  in-flight  refueling.  It  has  been  found  experimentally,  as  reviewed  in  [6]  and  [7],  that 
passive  surface  textures,  e.g.  streamwise  riblets,  can  reduce  the  turbulent  drag  on  a  surface 
by  five  to  ten  percent.  We  have  confirmed  this  numerically  and  have  shown  that  riblets 
work  by  damping  the  near- wall  spanwise  fluctuations  [8.9],  However,  we  are  unaware  of 
recent  efforts  to  control  or  constrain  the  lateral  spread  of  the  late  stages  of  natural  or  bypass 
transition  by  such  passive  texturing.  We  saw  a  unique  opportunity  here:  surface  textures 
like  riblets  have  not  seen  wide  application  due  to  their  relatively  modest  benefits  in  reducing 
turbulent  drag  (compared  to  their  weight,  installation  and  maintenance  costs,  etc.).  Locally, 
however,  where  the  flow  is  transitioning,  there  is  such  a  large  change  in  surface  shear  stress 
that  riblets  might  find  use  if  they  can  be  optimized  to  delay  or  constrain  such  transition. 

The  Force  Field  Model  and  Spectral  Approach:  There  exist  several  approaches  to  simu¬ 
lating  turbulent  boundary  layers  over  textured  surfaces.  We  have  published  a  method  for 
modeling  a  variety  of  surfaces  in  a  spectral  method  simulation  [8,9,10],  Our  approach  is  par¬ 
ticularly  flexible  and  has  been  successfully  tested  on  a  range  of  turbulent  flow  configurations 
over  riblets  and  MEMS-like  devices.  In  a  spectral  method,  the  spatial  distribution  of  the 
dependent  parameters  (velocity,  vorticity,  pressure,  etc.)  is  commonly  modeled  with  Fourier 
and  Chebychev  series  in  a  simple  domain.  Part  of  the  computation  is  performed  on  the  pa¬ 
rameters  themselves  and  part  on  the  Fourier  and  Chebychev  coefficients  of  the  parameters. 
Efficient  transform  methods  are  used  to  move  between  the  real  and  spectral  representations. 
The  spectral  code  we  use  is  based  on  the  method  in  [11]  (also  described  in  [12]). 

The  geometry  of  interest  here  is  that  of  a  developing  boundary  layer  in  which  the  wall  has 
an  array  of  riblets  or  fins.  The  spatial  development  in  the  normally  homogeneous  direction  is 
modeled,  as  in  Goldstein  et  ai  [13],  via  an  absorbing  buffer  layer,  Our  technique  for  modeling 
the  virtual  surface  introduces  a  localized  body  force  field  into  the  Navier  Stokes  equations. 


Such  an  approach  has  been  used  in  many  other  applications,  particularly  by  Peskin  [14]  and 
later  co-workers,  and  some  of  that  literature  was  briefly  summarized  in  [8]  and  [13].  Several 
groups  at  the  AFOSR.  LES-DNS  Summer  2001  meeting  discussed  related  approaches.  Our 
force  field  is  made  to  adapt  to  the  flow  and  bring  it  to  a  specified  velocity  on  the  intended 
boundary  points  thereby  creating  a  virtual  (stationary  or  moving)  solid  surface.  The  force 
is  computed  by  a  feedback  scheme  in  which  the  velocity  is  used  to  iteratively  determine  the 
desired  force  field.  The  stability  of  our  approach,  in  particular,  has  now  been  investigated  in 
detail  by  Lee  [15].  Our  force  field  technique  allows  us  to  utilize  the  fast  transform  methods 
without  encountering  the  severe  limitations  of  a  simple  geometric  domain  imposed  by  most 
other  spectral  approaches. 

This  virtual  surface  model  has  been  used  to  simulate  the  flow  around  moving  boundaries, 
vortex  shedding  off  of  circular  cylinders  [10]  and  cubic  obstacles,  the  turbulent  flow  over  flat 
and  riblet  covered  surfaces  [8,9],  arrays  of  wires  parallel  to  an  underlying  flat  surface  [8], 
arrays  of  linear  microjets  under  turbulent  boundary  layers  [16],  and  surfaces  having  three 
dimensional  waves.  Our  simulations  of  turbulent  flow  over  flat  plates  and  riblets  have  been 
validated  in  detail  against  both  experiments  and  other  workers’  simulations.  The  key  results 
from  the  riblet  work  related  to  the  present  effort  are  that  riblets  reduce  drag  by  damping 
the  cross-flow  fluctuations  near  the  riblet  crests,  that  some  riblets  cause  slugs  of  near- wall 
enstrophy  to  ride  along  the  rib  crests  parallel  to  the  mean  local  flow,  and  that  riblets  can 
turn  the  near-wall  flow  to  more  nearly  follow  the  riblet  crests. 

Proposal  Objectives:  Surface  Textures  to  Control  the  Growth  of  Spots 

I  proposed  to  continue  to  use  our  CFD  approach  to  investigate  the  use  of  surface  tex¬ 
tures  to  constrain  the  development  of  local  turbulent  regions.  As  a  turbulent  spot  moves 
downstream,  it  spreads  laterally  until  it  merges  with  other  growing  spots  and  the  entire  span 
becomes  turbulent.  Although  the  spots  themselves  have  an  interesting  structure,  they  appear 
much  like  small  arrowhead-shaped  regions  of  fairly  developed  turbulence  and  so  might,  be 
subject  to  means  of  control  already  developed  for  a  fully  turbulent  boundary  layer.  Although 
the  development  of  spots  may  be  unavoidable,  perhaps  their  growth  can  be  constrained  by 
a  riblet  covered  surface.  So  even  though  having  riblets  or  other  textures  covering  much  of 
an  aircraft  surface  does  not  usually  pay  off,  riblets  applied  locally  to  extend  the  end  of  the 
laminar  region  have  greater  potential  leverage  since  there  is  such  a  large  increase  in  shear 
stress  in  a  turbulent  flow  and  since  small  areas  of  riblets  are  cheap. 

Figure  1  indicates  a  schematic  idea  of  the  problem.  We  can  run  laminar  flow  over  a 
surface  covered  with  riblets  and  momentarily  introduce  a  localized  perturbation  to  induce 
turbulent  spots.  We  use  a  buffer  (fringe)  region  to  study  spatially  developing  flows  and  this 
requires  a  large  domain.  We  then  simulate  spot  growth  and  development  over  ribbed  and 
smooth  surfaces  to  investigate  the  differences.  We  have  examined  questions  such  as:  Do 
the  spots  grow  at  the  same  rates  over  different  surfaces?  Do  they  move  at  the  same  speed? 
How  are  the  internal  structures  of  spots  altered  by  the  surface  texture?  What  are  optimal 
riblet  shapes?  Do  riblets,  in  fact,  promote  the  growth  of  some  3D  instabilities?  (Presumably 
large  enough  riblets  do.)  Can  riblets  constrain  the  spreading  angle  of  the  spots?  Can  riblets 


Figure  1:  Schematic  showing  a  turbulent  spot  aft  of  a  perturbation  moving  over  riblets. 


usefully  change  the  recovery  to  laminar  flow  in  the  wake  of  the  spots? 

Our  initial  task  was  to  bring  the  new  student  up  to  speed  with  our  existing  software.  He 
was  then  to  examine,  separately,  baseline  cases  of  turbulent  and  laminar  flow  in  a  simple 
channel  and  over  riblets.  He  would  then  introduce  a  perturbation  in  a  laminar  Poiseuille 
flow  and  compare  to  previous  turbulent  spot  simulations.  Finally,  he  would  perturb  laminar 
flow  over  riblets.  The  results  are  presented  in  an  attached  detailed  report.  During  the  course 
of  the  project,  one  US  graduate  student  was  supported  as  well  as  two  US  undergraduates 
(both  of  whom  have  now  gone  on  to  graduate  school).  One  AIAA  conference  paper,  three 
technical  reports  and  one  MS  thesis  were  produced.  One  journal  publication  is  soon  to  be 
submitted.  The  graduate  student  spent  two  weeks  at  the  Navy  Research  Lab  in  Washington 
learning  about  the  spectral  code  and  has  now  completed  his  oral  and  written  PhD  qualifier 
exams.  Related  work  is  continuing  under  a  follow-on  AFOSR  grant. 

Detailed  Report: 

See  attached. 
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